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Computing Frequency by using Generalized Zero-Crossing Applied 
5 to Intrinsic Mode Functions 

ORIGIN OF INVENTION 

The inventor of the invention described herein is an employee of 
the United States Government. Therefore, the invention may be 
10 manufactured and used by or for the Government for governmental 
purposes without the payment of any royalties thereon or 
therefor , 

BACKGROUND OF THE INVENTION 

15 The need for instantaneous frequency (IF) is a real one for data 
from nonstationary and nonlinear processes. If the process is 
non-stationary, the frequency should be ever changing, albeit at 
a slow rate. Then, there is a need for frequency value as a 
function of time, for the value will not be constant throughout. 

20 For the nonlinear cases, the frequency is definitely modulating 
not only among different oscillations, but also within one 
period. A detailed explanation of IF is disclosed in related 
Patent Application Serial No. 10/615365, filed on July 8, 2003, 
entitled "Computing Instantaneous Frequency by Normalizing 

25 Hilbert Transform", inventor Norden E. Huang, which is 

incorporated by reference and assigned to the same assignee as 
this application. A computer implemented method of computing IF 
is disclosed in U.S. Patent Nos . 5,983,162, 6,381,559, 
6,311,13 0, all of which are also incorporated by reference. 
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The above-disclosed method of computing IF includes two 
essential steps and the associated presentation techniques of 
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the results. The first step is a computer implemented Empirical 
Mode Decomposition to extract a collection of Intrinsic Mode 
Functions (IMF) from nonlinear, nonstationary signals. The 
decomposition is based on the direct extraction of the energy 
5 associated with various intrinsic time scales in the signal. 
Expressed in the IMF's, they have well-behaved Hilbert 
Transforms from which instantaneous frequencies can be 
calculated. The second step is the Hilbert Transform of the IMF. 
The final result is the Hilbert Spectrum. Thus, the method can 

10 localize any event on the time as well as the frequency axis. 

The decomposition can also be viewed as an expansion of the data 
in terms of the IMF's. Then, these IMF's, based on and derived 
from the data, can serve as the basis of that expansion. The 
local energy and the instantaneous frequency derived from the 

15 IMF's through the Hilbert transform give a full energy- 
frequency- time distribution of the data, which is designated as 
the Hilbert Spectrum. 

However, there is a need for an alternative method to the 
20 Hilbert Transform in constructing an energy- frequency- time 

distribution. The Hilbert Transform is an intensive calculation 
process that requires powerful computers beyond many end users, 
including ordinary engineers and scientists, resources. 
Furthermore its results are not intuitive. If a signal to be 
25 analyzed has a surge-like behavior, i.e. the signal contains a 
very high amplitude for a short period time such as sound signal 
of a gunshot, the method disclosed in the above references may 
not generate reasonable results due to limitations imposed by 
spline fitting, which is used in the EMD in constructing extrema 
30 envelopes of the signal. The problem with spline fitting applied 
to the situation just mentioned, might result in both 
overshoot /undershoot problems and eventually divergence of 



IMF ' s . 



Zero Crossing points are the point where the voltage polarity 
of a waveform changes from negative to positive (or vice- 
5 versa) as it crosses the zero axis. In the field of signal 
processing, for a subject signal, the zero crossing points are 
detected and counted for analyzing the signal or making a 
decision for further actions. For example, zero crossing 
points are used to define bandwidth of a signal, which is 
10 assumed to be stationary and Gaussian. The bandwidth can be 
defined in terms of spectral moments as follows. The expected 
number of zero crossings per unit time is given by 

No = -{ — ) . 
^ mo 

where the expected number of extrema per unit time is given by 
1 ,m4s'^ 

15 A^. = i(^) ' 

^ 1712 

in which jfi. is the ith moment of the spectrum Therefore, the 
parameter, v, defined as 



20 offers a standard bandwidth measure. For a narrow band signal v 
= 0, the expected numbers of extrema and zero crossings have to 
equal. However, the previous equation defines the bandwidth in 
the global sense; it is overly restrictive and lacks precision 
at the same time as there might be more than two extrema between 

25 two zero crossing points. Consequently, the bandwidth 
limitation on the Hilbert transform to give a meaningful 
instantaneous frequency has never been firmly established. 
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SUMMARY OF THE INVENTION 

The present invention is directed to a computer implemented 
method of analyzing a signal. The method includes the steps of 
5 inputting the signal, extracting a set of Intrinsic Mode 
Functions from the signal, and generating a set of mean 
frequency functions from the Intrinsic Mode Functions. 

The present invention is directed to another computer 
10 implemented method of analyzing a signal. The method includes 

the steps of inputting the signal, extracting a set of Intrinsic 
Mode Functions from the signal, and generating instantaneous 
frequency based on critical points of the signal. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figures 1(a) -(b) are a high-level flowcharts describing the 
overall inventive method of extracting Intrinsic Mode Functions 
(IMFs) via Empirical Mode Decomposition (EMD) and using 
Generalized Zero Crossing to obtain the Instantaneous Frequency. 

Figure 2 is a graph illustrating an upper envelope, a lower 
envelope, an envelope mean and an original signal data. 

Figure 3(a) is a graph of a sound signal of a gunshot. 
25 

Figure 3 (b) is a graph describing the steps of constructing 
upper and lower envelopes of a part of the signal of Figure 3(a) 
with cubic spline fitting and generating an envelope mean from 
the envelopes . 

30 

Figure 3(c) is a graph describing the steps of constructing 
upper and lower envelopes of a part of the signal of Figure 3 (a) 
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with straight line fitting and generating an envelope mean from 
the envelopes. 

Figure 3 (d) is a graph displaying a collection of Intrinsic Mode 
5 Functions (IMFs) extracted via Empirical Mode Decomposition 
(EMD) with cubic spline fitting as in Figure 3(b). 

Figure 3(e) is a graph displaying a collection of Intrinsic Mode 
Functions (IMF) extracted via Empirical Mode Decomposition (EMD) 
10 with straight line fitting as in Figure 3(c). 

Figures 4{a)-{b) are graphs of data describing the Generalized 
Zero-Crossing approach. 

15 DETAILED DESCRIPTION OF PREFERRED EMBODIMENT OF THE INVENTION 

Figure 1(a) illustrates the overall inventive method of 
extracting IMFs via EMD including the Sifting Process in steps 
220 through 270. First, the signal can be inputted in step 
100 . 

20 

After inputting the signal in step 100, the analog signal is 
converted to the digital domain suitable for computer 
processing in the A/D conversion step 105. Depending upon 
whether the input signal is analog or digital step 105 may be 
25 bypassed. Thereafter, the Sifting Process (steps 107 through 
270) is applied to Sift the signal with the Empirical Mode 
Decomposition method and thereby extract the intrinsic mode 
function (s) . 

30 The signal from step 100 is first windowed by framing the end 
points in step 107. Then, the Sifting Process begins at step 
200 by identifying local maximum values of the digitized, 
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framed physical signal from step 107. 



Then, the method constructs an upper envelope 20 of the signal 
10 in step 210. The upper envelope 20 is shown in Figure 2 
5 using a dot-dashed line. This upper envelope 20 is 

constructed with a cubic spline that is fitted to the local 
maxima . 



Next the local minimum values of the signal 10 are identified 
10 in step 22 0. To complete the envelope, a lower envelope 3 0 is 
constructed from the local minimum values in step 230. The 
lower envelope 30 is also shown in Figure 2 using a dot-dash 
line. Like the upper envelope 20, the lower envelope 30 is 
constructed with a cubic spline that is fitted to the local 
15 minima. 

The upper and lower envelopes 20, 30 should encompass all the 
data within the physical signal 10. From the upper and lower 
envelopes 20, 30, an envelope mean 40 is the determined in 
20 step 240. The envelope mean 40 is the mean value of the upper 
and lower envelopes 20, 30. As can be seen in Figure 2, the 
envelope mean 40 bisects the physical signal 10 quite well. 



Referring to steps 300 through 330, Generalized Zero-Crossing 
25 (GZC) is a method of computing the Instantaneous Frequency for a 
nonlinear and nonstationary signal based on full and partial 
periods around a data point 8 of the signal as shown in Figure 
4. Zero crossing points are points at which a signal under 
consideration crosses a zero value. 

30 

In the GZC approach, all zero-crossings 2 and local extrema 4, 6 
are defined as critical points. As shown in the Figure 4(a), 
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the local extrema consist of local maxima 4 and local minima 6. 
For this approach, the signals to be analyzed are only meaningful 
Intrinsic Mode Functions (IMFs) when obtained via Empirical Mode 
Decomposition (EMD) . 

5 

The time intervals between all the combinations of the critical 
points 2, 4, 6 are considered as the whole or partial wave 
period [T^ , Tl , and 7^4). For example, the period between two 
consecutive up (or down) zero-crossings or two consecutive 

10 maxima 4 (or minima 6) can be counted as one period. Each given 
point along the time axis will have four different values from 
this class of period, designated as T^. Then, the period between 
consecutive zero-crossings (from up to the next down zero- 
crossing, or from down to the next up zero-crossing) , or 

15 consecutive extrema (from maximum 4 to the next minimum 6, or 
from minimum 6 to the next maximum 4) can be counted as a half 
period. Each point along the time axis will have two different 
values from this class of period, designated as T^. Finally, the 
period between one kind of extrema to the next zero-crossing 

20 point, or from one kind of zero-crossings to the next extremum 
can be counted as a quarter period. Each point along the time 
axis will have one value from this class of period designated as 
T4. Clearly, the quarter class, T^, is the most local and the 
full wave class, T^, is the least local. In total, there exist 

25 seven different period values at each point along the time axis. 
As T^,T2, and are weighted by factors of 4, 2 and 1, 
respectively, the weighted mean frequency at each point along 
the time axis can be computed as 



12 7; Tl Tl Tl 

wherein 
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CJ is mean frequency; 

are full periods (x = 1, 2, 3, and 4) enclosing the point 
under consideration; 

are half periods (y = 1 and 2) enclosing the point under 
consideration; and 

is a quarter period enclosing the point under 
consideration; 



An alternative is to compute a non-weighted as 



7 47; ITl 2T^ Tl Tj' 7]' T,' 

5 

It is to noted that various combinations of the critical points 
and weights can be used to compute the mean frequency. 



Referring to step 300, a mean frequency function 42 of the 
10 signal 41 can be obtained based on a collection of the mean 
frequencies along the time axis as shown in Figure 4 . This 
process continues until all of the mean frequency functions are 
obtained for all of the IMFs as in steps 3 00 through 3 05. 
Summing up the mean frequency functions (step 310) generates a 
15 function of frequency versus time, which is the instantaneous 
frequency. 

This approach is the most direct, local, and also the most 
accurate in the mean. Furthermore, this approach will also give 
20 a statistical measure of the scattering of the frequency value. 

This approach provides its crude localization, only down to a 
quarter wavelength. This approach tends to display its 
inability to represent the detailed waveform distortion; it 
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admits no harmonics. Unless the waveform contains asymmetries 
(either up and down, or left and right) , the GZC will give it 
the same frequency as a sinusoidal wave. For most practical 
applications, this mean frequency localized down to quarter a 
5 wave period is already a well-accepted result. As this method 
physically measures the period, or part thereof, the values 
obtained can serve as the best local mean over the period to 
which it applies. 



10 Referring to step 350 of Figure 1(a), serious problems of the 
spline fitting can occur near the ends, where the cubic spline 
being fitted can have large swings. The end swings, if not 
fixed, can eventually propagate inward and corrupt the whole 
data span especially in the low frequency components. For 

15 example, if the signal contains a very high amplitude for a 

short period time such as sound signal of a gunshot as shown in 
Figure 3 (a) , the method disclosed in the above references many 
not generate reasonable results due to extreme overshoot in the 
cubic spline fitting. Referring to Figure 3(c), which deals 

20 with a problematic area of the data of Figure 3 (a) , the upper 
and lower envelopes 21, 31 are constructed by connecting local 
maxima and minima with straight lines, respectively. For 
comparison. Figure 3(b) shows the upper and lower envelopes 21', 
31', which are constructed with the cubic spline fitting. As 

25 anticipated, the EMD based on the cubic spline fitting does not 
generate reasonable IMFs as shown in Figure 3 (d) whereas IMFs 
obtained through the straight line fitting gradually converge to 
a monotonically decreasing signal as shown in Figure 3(e). 



30 It is to be understood that the above- described arrangements are 
only illustrative of the application of the principles of the 
present invention. Numerous modifications and alternative 
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arrangements may be devised by those skilled in the art without 
departing from the spirit and scope of the present invention and 
the appended claims are intended to cover such modifications and 
arrangements . 
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